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Abstract 

Two different methods of tlie covariant relativistic separable kernel construction 
in the Bethe-Salpeter approach are considered. One of them leads in the center- 
00 , of-mass system of two particles to the quasipotential equation. The constructed 

^ I 4-dimensional covariant functions are used to reproduce neutron-proton phase shifts 

CN I for total angular momenta J = 0, 1. Obtained results are compared with other 

^ ■ model calculations. 

^ , 

O ■ 1 Introduction 
m 

The relativistic description of the interaction between two nucleons can be carried out us- 
ing the Bethe-Salpeter (BS) equation [Ij and different quasipotential approximations to 
it [2]. One of the ways to solve the BS equation is using the separable presentation of 
its interaction kernel. The BS formalism with such kernel have been successfully applied 
^ i to a covariant description of the elastic lepton-deuteron scattering, deuteron electro- and 
photodisintegration, deep inelastic scattering (DIS) of leptons on light nuclei. It has also 
provided good results for electromagnetic reactions with the deuteron. The approach facili- 
tates the analysis of the deuteron properties and its comparison with nonrelativistic results 
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^ . (see, review [3] and references therein). 
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The consideration of the deuteron break up reaction demands the final state interaction 
between the outgoing nucleons to be taken into account. If the separable kernel of interac- 
ts ■ tion is used it leads to the calculation of integrals with separable functions. Usually they 
^ ■ represent a generalization of nonrelativistic separable functions [1]. The presence of second- 
and higher-order poles in the integrand containing these functions makes such calculation 
impossible. In the paper [5] we constructed the new parametrization with functions sug- 
^ I gested in [6|, which do not contain poles on the real axis in the relative energy complex 
d ' plane, and apply them to the description of the neutron-proton (np) scattering. Within 
this approach we have achieved a good agreement with experimental data in the region of 
high kinematic energies of the nucleons TLab ~ 1-2 GeV. But as it was resumed in ^ it is 
necessary to improve an agreement with experimental data in a whole energy region up to 
^Lab < 3 GeV by increasing the rank of the separable kernel of the nucleon-nucleon (NN) 
interaction. This work is in progress. 

Before that we investigate a case when the separable NN kernel which represents a 
3-dimensional function in the center-of-mass system (CMS). The use of such functions 
leads to the approaches called quasipotential where the two-body Green's function in the 
BS equation is replaced by a suitable function which allows to reduce the 4-dimensional 
equation to 3-dimensional one, for example, like in |2]. 

In this work the one-rank separable kernels of NN interaction for the partial states with 
total momenta J = 0, 1 are presented for functions in 4-dimensional and 3-dimensional 
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forms. Parameters of the kernels are fitted by the calculation of phase shifts and low- 
energy observables (scattering length, effective range, deuteron bound state energy) using 
the constructed T matrix. 



2 Methods of a covariant relativistic generalization 

One of the common methods of a covariant relativistic generalization of the Yamaguchi- 
and Tabakin-type functions is to replace 3-momentum squared by 4-momentum squared: 

2 2 2,2 /i\ 

P ^ -P =-Po + P ■ (1) 

This formal procedure converts 3-dimensional covariant functions to 4-dimensional ones. 

Another method is based on the introduction of the formal 4-vector Q though the 
relative p and total P 4-momenta of the two-body system by the following relation: 

Q=P-^P, (2) 

with total momentum squared s = P^. 

Note that in CMS where P = {^/s, 0) the 4-vector Q is defined by components Q = (0, p) 
and thus 

P' = (3) 

can be formally converted to the Lorentz invariant. 

Let us consider the methods described above in application to the nonrelativistic Yamaguchi- 
type function 

9i\p\) = (4) 
In the first case using the substitution ([1]) we obtain the covariant function in the form: 

In the second case we use the relation ([3]) and obtain the function: 

The presented functions have rather different properties in the relative energy po complex 
plane in CMS. The function Qp has two poles on the real axis for pq at ±Vp^ + P'^ -F 
while the function qq has no poles on it. 

In practical calculations of the reactions with the high momentum transfer the po inte- 
gration can lead to singular expressions in functions Qp^Q on |p| or cos 6p. Such problem 
can be easily solved by calculating the \p\ or cos 6p principal value integral. However, 
another form of functions Qp^q with odd powers in a denominator leads to nonintegrable 
singularities. Therefore we introduce functions Qp^Q of type without poles on the real axis 
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of the relative energy pq complex plane. As an example of such function we introduce the 
covariant form factors in the following form (see also the section 3 of the paper [5]): 

n(vP) = ^ Pc-pI+P' .7. 

and in the second case we use the relation ([3]) and obtain the function 

a (n P)= m P^+P' (H) 

UQKP, ) _ + ^4 (p2 + ^2)2 + ^4 • ^ ^ 

We note that the function qq still has no poles on the Pq real axis while the Qp has poles at 
Po: ± Vp' + + ±Vp' + f3'' - ia'. 

The two methods of a covariant relativistic generalization described above can be inves- 
tigated by solving the Bethe-Salpeter equation for specific partial states. 



3 Solution of the Bethe-Salpeter equation 

The details of the solution of the BS equation with the one-rank separable kernel of inter- 
action can be found in the paper [5]. The value to be calculated is h{s): 

his) = J dpo J \p\'^d\p\S{po, \p\; s)g{po, \p\f (9) 

with 

being the two-particle Green's function; here Ep = y/p'^ + m'^, m is the nucleon mass. 

To obtain the function h{s) the 2-dimensional integral onpo and \p\ should be calculated. 
To perform the integration on pq the Cauchy theorem is used. As it can be seen from 
Eq.dnD there are two types of singularities in the Pq complex plane: first ones are poles of 
the function S{po, \p\; s): 

p^^'^^ = ±yrs/2TEpTte, (11) 

and other - poles of the function g{po, \p\). 
The function Qp has four poles: 



= ±Jp^ + ia\ (12) 



and to perform the po integration residues in three poles of Eqs. ljTTjl and (1121) should be 
calculated. These calculations are performed analytically. 

The function gq has no poles on the po real axis and therefore the only poles of Eq. ffTTl) 
should be taken into account. The result for h{s) can be written as: 

Ks) = \p\Mp\ • (13) 

^ ' 27r2 i 7i - 2Ep + ie ^ ^ 
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This equation formally coincides with that could be obtained within the Blankenbeckeler- 
Sugar-Logunov-Tavkhelidze (BSLT) approximation [2] which amounts to replacing the 
Green's function in Eq.([9]) by the expression 

SbslApo, \p\; s) = ~2m{^ - 2Ep + te)-^S{po). (14) 

Although the solutions of the equation with functions qq and in the BSLT approximation 
coincide in CMS the difference becomes evident when the reaction with the two-particle 
system is considered. In that case the arguments of the function qq are calculated with the 
help of the Lorentz transformations in the system different from CMS. 



4 One-rank kernel 

We analyze two covariant relativistic generalizations of the Yamaguchi form factors: mod- 
ified Yamaguchi (MY) functions and modified extended Yamaguchi (MEY) functions. 

4.1 Modified Yamaguchi functions 

For the description of the chosen partial states we use the following covariant expressions 
for both cases: 

9p IP0,P)- ^^2_p2_^2)2 + ^4' (15) 



_j_ 



[S], N [51/^ CMS (Pcl+P^) /.-s 



91%0,P) = .r(Qo, Q) ^ (18) 

It should be noted that the 4- vector Q is defined by the equation ([2]). The functions are 
numerated by angular momenta L = 0{[S]),1{[P]). The numerator in g^q is introduced to 
compensate an additional dimension in the denominator to provide the total dimension as 
GeV~^. This form was chosen because at pci = Z^^, « = we get the function g^^ g in the 
standard Yamaguchi form |7j. We prefer not to consider the case with the square root in 
the denominator as it is used in [S] for S waves to avoid the calculations with cuts on the 
real axis in the po complex plane. 
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Table 1: Parameters of the one-rank kernels with MYQ functions. 



T^max 
-'Lab 


1 Q + 
Oq 

0.999 GeV 


■'St 

1.1 GeV 






1.1 GeV 


3 p+ 

0.999 GeV 


3p+ 

0.999 GeV 


A, GeV2 


-351.73 


-0.048442 


A, 


GeV^ 


9.2481 


-172.92 


14.245 


f3, GeV 


1.8805 


0.19973 




GeV 


0.56446 


0.72278 


0.58444 


a, GeV 


0.24919 


1.1891 


a, 


GeV 


0.2 


0.68577 


0.2 


Pel, GeV2 


-4.6069 


-332.95 












So, GeV2 


4.0279 


4.2020 


3.8682 



Table 2: Parameters of the one-rank kernels with MEYQ functions. 



T-imax 
-'Lab 


1 C+ 

0.999 GeV 


''St 
1.1 GeV 






1.1 GeV 


3 p+ 

2.83 GeV 


3p+ 

0.999 GeV 


A, GeV^ 


-0.95380 


-16.852 


A, 


GeV^ 


0.02 


-197.47 


0.031096 


Ci2, GeV" 


3.2264 


18.606 


Ci2, GeV 


0.15546 


6.2551 


-0.071656 


(3u GeV 


0.10621 


2.0098 


A, 


GeV 


0.16353 


0.67649 


0.19350 


p2, GeV 


0.1 


0.1 




GeV 


0.44029 


1.2192 


0.65669 


ai, GeV 


1.4281 


0.20562 


ai, 


GeV 


0.2 


0.74382 


0.2 


a2, GeV 


1.9491 


1.1581 


"2, 


GeV 


0.53818 


1.1973 


0.2 


Pel, GeV2 


27.424 


132.22 


Pe3 


, GeV2 


103.92 


-2.8247 


-630.35 


Pe2, GeV2 


10.690 


-0.293 













4.2 Modified extended Yamaguchi functions 

To extend the form of functions with increasing number of parameters we introduce the 
following gp Q functions: 

[S]. . _ (Pel -pI+ P") Cl2(pg - P'){j>e2 -pI+ P"? , . 

{pi - - + «f {{pi -p'- Ply + c^iy ' ^ ^ 



„) = V -^°+^^ + c,2^i-Pi+P^nPcs-Pi+P^) 

[^1^ ^ [S]fn n\ (Pci + p"^) Ci2P^{p^2+P^? (^^. 

{Po,P) = i KQo,Q) (p2 + ^2)2 + ,4 - (pqr^|)iT^' (21) 

[P]( X [P]/^ ^^ CMS IpI , Ci2\p\^{pe3+P^) 

The denominators with Pci,Pe2,Pc3 are intended for a dimension compensation as Pd for 
g^pj^ in the previous case (fT^ . (fT7|) . 
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Table 3: The low-energy scattering parameters and binding energy for the singlet (s) and 
triplet {t) S waves. 





as(Fm) 


ros(Fm) 


at(Fm) 


rot(Fm) 


^d(MeV) 


MYQ 


-23.69 


2.22 


5.37 


1.71 


2.2246 


MEYQ 


-23.49 


2.30 


5.37 


1.71 


2.2246 


Experiment 


-23.748(10) 


2.75(5) 


5.424(4) 


1.759(5) 


2.224644(46) 
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Figure 1: Phase shifts for the ^P^' wave. Dash (short dash) line corresponds to our 
parametrization with MY (MYQ) functions, eq. (fT6!) ( (fT8!) ). respectively; the solid (dash 
dot dot) line illustrates the extended case MEY (MEYQ), eq.([20])(([22D); the dotted and 
dash-dotted lines describe the CD-Bonn and SP07 calculations, correspondingly. 

5 Discussion 

Using the np scattering data we analyze the parameters of the separable kernels as it was 
earlier presented in [5]. 

The calculated parameters of the one-rank kernels with MYQ and MEYQ functions 
and Tj™^^ are listed in Tables 1 and 2. In Table 3 the calculated low-energy scattering 
parameters for S waves are compared with experimental data. 

The calculated phase shifts are presented in FigsHHS] where MY (MYQ) case is denoted 
by the dash (short dash) line, MEY (MEYQ) one - by the solid (dash dot dot) line. Parame- 
ters for the first case (functions MY and MEY) have been presented in [5] . The calculations 
with the nonrelativistic CD-Bonn [8] (dotted line) potential and the SP07 solution [9] (dash 
dot line) are included for comparison. 
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In FigHlwe can see that all the calculations give a reasonable agreement with experi- 
mental data up to the energy TLab~ 1-1 GeV but at larger energies their behavior becomes 
different drastically. Note that the CD-Bonn and SP07 results for phase shifts even change 
sign at 1.5 < TLab < 2.5 GeV. So we can stress that it is very desirable to obtain more exact 
determination of experimental data for the ^P^ channel in this energy range. Now keeping 
in mind that the discrepancies between CD-Bonn, MEY, MEYQ and SPOT calculations at 
^Lab < 1.2 GeV are not so large we can trust our calculations (in both cases) in this range 
of energies. 




Figure 2: The same as in Fig.l for the ^Pg^ state. 



In Figl2] CD-Bonn and MY (MYQ) calculations demonstrate an opposite behavior for 
^Pq' phase shifts at TLab > 1 GeV. The MEY (MEYQ results practically coincide with 
MEY) and SP07 calculations give a very good agreement with experimental data in a wide 
range of energies TLab < 3 GeV. Thus, both variants of the kernels of the NN interaction in 
the ^Pq~ channel are acceptable to be used in various relativistic calculations of reactions. 

Calculations for the ^Pi' channel (see FigjS]) show us that MY (MYQ), CD-Bonn, and 
MEY have a similar behavior in a wide energy range, but are able to explain experimental 
data up to TLab =0.55, 0.6 and 1.2 GeV, respectively. It should be noted that MEYQ 
calculations describe phase shifts up to 1.5 GeV. The SP07 calculations give a suitable 
agreement with experimental data up to 3 GeV. Calculations for the cases 1 (MY, MEY) 
and 2 (MYQ, MEYQ) coincide up to 0.55 GeV, but at larger energies their behavior becomes 
significantly different. So it is clear that to explain the data for the ^Pi' channel in a wider 
energy range within our approach we need to increase the rank of the separable kernel of 
the NN interaction. 
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T, ^, GeV 

Lab' 



Figure 3: The same as in Fig{T]for the ^P^^ state. 



FiglH demonstrates a close agreement with experimental data for CD-Bonn up to the 
energy TLab < 0.6 GeV, MY and MEYQ - up to TLab < 1.3 GeV in the ^S^ channel (note 
that these calculations give a very similar behavior of phase shifts in a wide energy range 
up to 3 GeV). The MEY calculations give a very good description of phase shifts up to 
1 GeV. So we can conclude that in this case as well as for "^P^ we need to increase the rank 
of the kernel. The SP07 calculations describe experimental data for the -^Sq channel up to 
3 GeV. 

FigjS] for the ^Si wave shows us that CD-Bonn calculations agree with experimental 
data up to 0.6 GeV. All the other compared results demonstrate a similar behavior in a 
wide energy range and explain known experimental data up to TLab = 1.1 GeV. It should 
be noted that calculations of low-energy parameters with MY (MYQ) and MEY (MEYQ) 
form factors give us a reasonable agreement with experimental data (see [5j for the first 
case and Table 3 - for the second). 



6 Conclusion 

We present two new different parameterizations (4-dimensional and 3-dimensional in CMS) 
of the separable kernel of the NN interaction which are adopted for calculations at large 
energies. As it was expected the MEY(MEYQ) functions give a better description of the 
scattering data than the MY(MYQ) ones. Using the suggested MEY (MEYQ) form factors 
the phase shifts are described in a whole range of measured energies for the following partial 
states: ^P^, ^St (Ttab < 1.2 GeV) and 3p+ (TLab < 3 GeV). The phase shifts for the ^S^ 
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I 1 I 1 I 1 1 
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Figure 4: The same as in Fig{T] for the singlet partial state ^5*0^, but the description by 
MY (MYQ) functions is defined by eq.(^{(^) and extended one MEY(MEYQ) - by 
eq.dlKdni)). 



and ^Pi' partial states can be described in our approach up to TLab < 1-1 GeV. To improve 
their agreement with experimental data up to TLab < 3 GeV it is necessary to increase the 
rank of the separable kernel of the NN interaction. This work is in progress. Finally we can 
say that 4-dimensional (Bethe-Salpeter type) as well as 3-dimensional (quasipotential type) 
in CMS parameterizations give a good description of the considered shift phases. Both of 
them are to be used in calculations of various electromagnetic reactions. 
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Figure 5: The same as in FigjUfor the triplet partial state ^Si. 
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